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ABSTRACT 


Detailed hydrodynamic modeling of the passage of supernova shocks through 
the hydrogen envelopes of blue and red progenitor stars has been carried out to 
explore the sensitivity to model conditions of light element production (specif- 
ically Lx and B ) which was noted by Dearborn, Schramm, Steigman and 
ruran (1989) (DSST). It is found that, for stellar models with M & 100M©, 
current state of the art supernova shocks do not produce significant light element 
yields by hydrodynamic processes alone. The dependence of this conclusion on 
stellar models and on shock strengths is explored. Preliminary implications for 
Galactic evolution of lithium are discussed, and it is suspected that intermediate 
mass red giant stars may be the most consistent production site for lithium. 


Introduction 

In recent years, lithium has become a very important probe of cosmological 
and Galactic history. As we will see, lithium yields are very sensitive to nucle- 
osynthetic conditions. Small changes in the environmental parameters can lead to 
large changes in the amount of lithium produced or destroyed by these processes. 
Because of this, it has become common to try to enhance our understanding of 
stellar and early Universe physics by using current lithium abundances to set lim- 
its on such astrophysical parameters as the average mass density in baryons of the 
Universe (Kawano, Schramm, and Steigman 1988) and the depth of convection 
zones in proto-stars (D’Antona and Mazzitelli 1984). A big uncertainty in this 
approach is that there are at least two distinct abundance measurements for 7 Li 
. The abundance of 7 Li in low metallicity, Population II main sequence stars in 
the Galactic halo is constant above a certain effective stellar temperature, with 
an observed (Spite and Spite (1982)) 7 Li /H ratio of KT 10 . Since the metals 
are assumed to have been produced as the Galaxy evolved and developed a disk, 
this group of stars is generally believed to be the earliest group of observed stars 
to have formed . This essentially constant value of the lithium abundance, over 
a small (but significant) range of stellar masses (all but the least massive pop II 
stars have long since left the main sequence) and a wide range of met alii cities, 
has reinforced the conclusion that this lithium abundance is primordial. This 
is consistent with standard big bang nucleosynthesis (BBN). However, the high 
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metallicity (z > .01) stars of the Galactic disk have a consistent flat upper limit 
(Boesgaard 1990) of 7 Li /H & 10~ 9 . This maximum abundance remains high 
even for the oldest pop I clusters (e.g. NGC188) (Hobbs and Pilachowski 1988). 
Thus it appears that 7 Li abundances grew along with the other heavier elements 
to some extent. 

In general, any changes in the early Universe scenario can only make the 
predicted primordial lithium abundance go up. This is because lithium is formed 
by two different nuclear processes: either as beryllium at high densities (via 
3 H e(a , 7 ) 7 Be) which then decays to lithium by electron capture with a lifetime of 
about forty-nine days, or as lithium directly (via ‘ i He(a,p)‘ Li) at low densities. 
The reason that 7 Li can only be made directly at low densities is that it is 
very easily destroyed by the 7 Li(p,a) A He reaction at what are relatively low 
temperatures (~ 2 x 10 6 °K ) due to the extreme stability of 4 He . (At stellar 
densities, lithium is rapidly destroyed, at temperatures above 10 6 °K . ‘Be on 
the other hand, cannot be destroyed by protons at temperatures high enough to 
create it, due to photo-disintegration reactions (see below)). The gap between 
these two cosmological production regimes causes a dip in the lithium abundance 
versus density of the Universe curve (figure 1). The standard BBN abundance 
brackets the bottom of the dip at Li 7 / H ss 10 -1 °. All other theories of primordial 
nucleosynthesis (for example inhomogeneous nucleosynthesis (Applegate, Hogan, 
and Scherrer 1987) or late decaying particle nucleosynthesis (Dimopoulus, et. al. 
1988) predict higher primordial values of lithium than this, since any change in 
conditions will favor one or the other production mechanism to a greater degree 
than the standard scenario and this will cause a greater primordial 7 Li yield. 

If the pop II abundance is primordial, then we need to make ~ 10 times 
our starting abundance in some Galactic process by the time of the formation 
of the earliest pop I clusters. Since making lithium directly at stellar densities 
is impossible (any temperature hot enough to make ‘Li will immediately de- 
stroy it at a higher rate), any stellar lithium production process must follow the 
Cameron-Fowler (1971) paradigm of making 7 Be in a hot environment and then 
cooling it before the ~ 49 day half-life runs out. The original Cameron-Fowler 
proposal used a convective red giant envelope as the production site. However, 
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DSST noted that similar conditions might also be achieved in supernovae. Since 
supernovae are by far the most theoretically favored sources of the metals, it 
seems reasonable that they might also be the sources of the 7 Li . This would 
explain why the 'Li abundance “tracks” the metal abundance. The traditional 
Cameron- Fowler (1971) site of red giants and/or the alternative of classical novae 
(Starrfield, Truran, Sparks, and Arnould (1988)) might have difficulties explain- 
ing the pop I 7 Li abundance observed in the oldest pop I clusters because these 
sites involve lower mass stars which have longer lifetimes, and therefore might not 
put their produced lithium back into the interstellar medium rapidly enough. An 
alternative supernova source for 7 Li is the v-process of Woosley, et. al. (1990). 
However, there is some concern that the ^-process may tend to overproduce 11 B 
as well as rare odd-odd nuclei such as 180 Ta if it yields significant 7 Li. 

A related question to all of this is that of the n B abundance. BBN models 
do not make 11 B . It is generally assumed that 10 B and 6 Li are made by Galactic 
cosmic ray spallation of heavier nuclei (Reeves, Fowler, and Hoyle (1970)). How- 
ever, cosmic ray spallation produces (Walker, Matthews, and Viola (1985)) 7 Li 
and 6 Li in the ratio 7 Li/^Li ~ 2 whereas solar system observations (Cameron, 
1982) yield 7 Li/^L i ~ 12 in agreement with the lower bound of ' Li/ 6 Li ^ 10 
in F and G stars determined by Anderson, Gustafsson and Lambert (1984); for 
the interstellar medium towards Zeta Ophiuchus, Ferlet and Dennefeld (1984) 
find ‘Li/^Li ~ 38 (<J 25). Spallation also produces U B in an isotope ratio 
n J?/ 10 f? ~ 2.5 whereas the observed (Cameron 1982) boron isotope ratio is 

U B/ 1{) B ~ 4. Thus, if spallation produces 6 Li and 10 B (for which there are no 
other proposed mechanisms, since both 6 Li and 10 B seem too fragile for stellar 
production), it is unlikely that spallation can account for all 7 Li and 11 B produc- 
tion as well. Since the limiting reaction for 7 Be production in Cameron-Fowler 
schemes is 7 Be (a, 7 ) n C — > e”+ n i? we might suspect that the supplemental 7 Li 
and 11 B sources could be identical. 

In 1973, it was proposed by Colgate (1973) and Hoyle and Fowler (1973) that 
deuterium and other light elements might be created if a strong ion shock were 
produced in the very low density (~ 10“ 5 ^/cm 3 ) region of the hydrogen envelope 
of a core collapse supernova. Epstein, Arnett, and Schramm (1976) noted that 
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in this case, 7 Li and 11 B were particularly easy to make in the ion shock process. 
Since lithium and boron are rare, this limited the process, if it occurred at all, to 
the production of only 7 Li and 11 B and ruled out D production in this process. 
Weaver (1976) showed that for envelope shock models of up to very high energies: 


for E s hock < 100 MeV /nucleon = mfjVg hoc */2 -+ v 3 h 0C k = 1-4 x 10 10 cm/s, (1) 

(where m H is the mass of a hydrogen atom) the shocked gas remains essentially 
in radiative equilibrium and thus such shocks are mediated almost entirely by ra- 
diation (i.e. there was no significant ion shock). This implied that in the relatively 
diffuse medium of a red giant envelope, the shock-induced temperature jump was 
too small by the time the shock reached the base of the surface convection zone, 
where 3 He could be found to use as a starting material for 7 Be synthesis, to 
produce significant 7 Li . However, Amould and Norgaaxd (1975) showed that 
for certain generic explosive nucleosynthesis environmental conditions (we will 
discuss this scenario below), significant 7 Li production could occur. With the 
explosion of SN1987A, attention turned to the specific dynamics of core collapse 
in blue giant stars. DSST argued that in these stars, with their more compact 
envelopes, shock temperatures sufficient for 7 Li and n B production might oc- 
cur around the base of the 3 He zone, that is, processing conditions utilized by 
Arnould and Norgaard might occur (see figure 2). However, they noted the sensi- 
tivity of their conclusions to the shock assumptions. They also admitted that the 
shock physics used in their calculations might be less detailed than necessary for 
accurate conclusions. In this paper, we have used a state of the art calculation of 
the envelope shock dynamics to explore the sensitivities of the models. In partic- 
ular, we note that our current, more detailed treatment of the shock puts more 
energy into radiation pressure than the model of Dearborn et. al. This results 
in a lower temperature behind the shock front for a given configuration and thus 
less light element production. Woosley, et. al. (1990) have not only proposed 
a different mechanism (i.e. neutrino induced reactions) which could produce 7 Li 
and 11 B deeper in the supernova than the zones we consider but also questioned 
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whether the results of Dearborn et. al. would be confirmed by detailed model- 
ing, since the post neutrino shock wave processing calculations in their models 
produced significantly smaller light element yields. 

In the rest of this paper, we first discuss the models used to create the pre- 
supernova stars. We then give a detailed description of each of the components of 
the simulation of the envelope shock: the underlying Von Neumann- Richtmeyer 
artificial viscosity hydrodynamics code, our method of treating radiation pres- 
sure and transport, the approximations we assumed in modeling the fundamental 
physics of the gas of the stellar envelope, the “piston” which we used to drive the 
explosion, and our thermonuclear reaction network. Next, we discuss our meth- 
ods for checking the reliability of the code. Finally, we present our conclusions 
about the supernova shock yields of 1 Li and 11 B and discuss their significance 
in the light of the overall picture of Galactic chemical evolution and of other 
proposed production mechanisms. 


The Input Models 

Figure 2 shows a schematic of a typical stellar model immediately before it 
explodes. There is a very dense (~ 10 4 g/cm 3 ) carbon/oxygen core surrounded 
by a carbon shell (~ I0 3 g/cm 3 ), a helium shell (~ 100 — I0g/cm 3 ), and an 
extended hydrogen envelope containing 60 — 65% hydrogen by mass. There is a 
surface convection zone which mixes the material in the outer envelope down to 
the- region where the H and 4 ife mass fractions are approximately equal. This 
convection brings original 3 He (he. the 3 He that the star was “born” with plus 
the primordial 2 H that had burned to 3 H e during its lifetime) down to this point. 
The increased temperature will destroy most of the 3 He present below the base of 
the convection zone, but a substantial amount ( X « 10 -4 ) of 3 He will be present 
at this point. As the supernova shock weakens in its passage through the outer 
envelope, the temperature immediately behind the shock (the region where our 
proposed synthesis will occur) decreases such that the closer the innermost 3 He 
containing zone is to the center of the star, the more ‘Li should be produced. We 
expect that blue progenitor stars, with their more compact outer envelopes, will 
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produce considerably more 7 Li than red stars. The higher densities in the blue 
stars which lead to generally higher post-shock temperatures should also favor 
them as 7 Li production sites. These conjectures are bom out by our simulations. 

Specifically, the supernova precursor models we used came from the evo- 
lution of 15, 25, and 100 solar mass stellar models which started with initial 
compositions of a H mass fraction of .7, a metallicity z = 0.02 (standard Pop 
I metallicity), and a z He mass fraction X( z He) = 2.1 x 10~ 4 . This z He value 
assumes that all initial deuterium is processed to z He during the protostellar col- 
lapse phase (Dearborn, Schramm, and Steigman 1985). The models were evolved 
from the zero age main sequence into the core carbon burning phase, after which 
no additional nucleosynthesis occurs in the envelope prior to the core collapse. 
Since we will only be concerned with processing in the outer envelope, the mod- 
els were not evolved past this point. A more complete description of the stellar 
evolution code used can be found in the appendix to Dearborn, et. al. (1990). 

The radii of red giant models are sensitive to the surface opacities and the 
mixing length. For this reason, the mixing length is often taken as a free pa- 
rameter to fit the observed temperature of stars in a cluster. In our red giants 
(models ml50 and ml51), we used models with slightly different ages, as well 
as different mixing lengths (Table 1). While the radius of a star at the end of 
its lifetime is quite uncertain, we can be guided by comparison of our models 
with observations. Brown, Tomkin, and Lambert (1983) examined the available 
observations of Betelgeuse and found it to be a 15 solar mass star with a radius 
between 4.26 and 4.67 x 10 13 cm. Since it is a variable, the precise radius at the 
time of its core collapse cannot be pre-determined. Our models ml51 and ml52 
are, however, quite representative of Betelguese. The envelope contains 7.2 solar 
masses of hydrogen, and 11 solar masses total. The models are convective in the 
outer 10.5 Mq, and the z He mass fraction is a constant through this region, at 
approximately 10~ 4 . Below the base of the convection zone, the abundance of 
z He falls discontinuously to its equilibrium value outside the Hydrogen burning 
shell. 

Although our stellar models go through a “blue loop” period during Helium 
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burning, they all evolve to the red before the core collapse occurs. It is known 
that the outer hydrogen envelopes of giant stars can alternate between extended 
(red) and compact (blue) phases regardless of the evolutionary state of their 
cores, so that they may experience core collapse in either red or blue states. In 
particular, SN1987A exploded from a blue progenitor. In order to obtain a “blue” 
envelope structure that is similar to that which SN1987A must have had, we took 
models at a slightly earlier evolutionary phase, the end of helium burning and 
the beginning of carbon burning. Our models at this stage were still on a “blue 
loop” and had not yet returned to the red. Lithium production in such a star 
requires it to have had a red giant stage with its attendant convection zones, in 
order to distribute the surviving 3 He to the deeper regions of the outer envelope. 
Following this, the blue loop causes some additional 3 He depletion in the deeper 
regions. 

In order to cover a wide range of possible input models, we also used a model 
of a 100M© star in two configurations. We used one model with absolutely no 
mass loss (almost certainly an incorrect model), and one with mass loss included 
(final pre-explosion weight 58M©). Since our experience indicated that we could 
expect much greater 7 Li production in bluer models, and since we wanted to 
explore whether these extremely massive stars would show substantially different 
behavior from the lower mass models, we exploded these two models while they 
were in their “blue loop” state. 

Since 3 He is not present below the outer edge of the 4 #e zone, this was 
a convenient point at which to initiate our simulation. The outer envelope, 
starting at the outer edge of the A He zone, was “clipped off” and a “piston” was 
attached to the innermost zone. The “piston” driver which we used to simulate 
the behavior of the ejected material coming from the core of the star was modeled 
by fixing the velocity of the innermost zone to a specific predetermined function 
of time (see below). For modeling of the piston velocity, it is necessary to describe 
the outer hydrogen envelope as a set of regions of power law density decrease such 
that 


8 



— LU 


( 2 ) 


where p is the density, a a constant, r the radius, and w the powerlaw exponent. 
We chose regions by inspection from a log-log graph of p vs r (see figure 3a ). As 
can be seen, the envelope can be conveniently subdivided into a small number 
of regions of nearly exact powerlaw behavior. We found omega by a chi-squared 
linear regression fit for each “hand picked” region. Our linear correlation coeffi- 
cients ranged between .97 and 1.0, where 1.0 is a “perfect” line. Alternatively, 
our standard (see, for example, Press, et. al. (1988)) estimated errors m slope 

(u j) ranged from 0.3 — 5%. 

The zoning structure used for the stellar evolution code is not very useful 
for our purposes. The zones axe much too thick for the detailed examination of 
the primary area of interest centered around the base of the surface convection 
zone, we rezoned the models using the following prescription. All zones of mass 
10 30 - 10 32 g were not changed. Zones of mass > 10 32 <7 were cut into an integral 
number of equal mass zones, such that their masses were as close to I0 32 g as 
possible. Zones of mass < 10 3 % were grouped together into a single zone until 
that zone’s mass was > 10 30 </. Since this new zoning was not in strict hydrostatic 
equilibrium and since it resulted in adjacent zones of up to 50% difference in 
mass, in higher energy shocks there was some slight (£ 10%) “ringing of the 
hydrodynamic variables behind the main shock due to small shocks traveling 
between the real shock front and the beginning of the rezoned region. 


The Simulation 


Hydrodynamics 

The hydrodynamic code we used was a modified version of the SIMPLE code 
(Larsen 1989). This code is a one- dimensional (spherically symmetric, cylindri- 
cal^ symmetric, or slab geometries can be selected) Lagrangian mesh model. 
The SIMPLE code treats the gas in which the shock takes place as three sepa- 
rate fluids (electrons, ions and radiation) and links them together with simplified 
models of the relevant physics. We ran the simulations in “full ionization” mode, 
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t.e. assuming that the nuclei were completely ionized at all times (essentially 
exact, since the lowest temperatures we dealt with were in excess of 10 6 °I\). 
Since we did not run with full radiation transport (see below), this ionization 
assumption did not affect the opacity of the stellar material but only affected the 
hydrodynamics and heat diffusion parts of the code. The electrons, the ions, and 
the radiation were thus considered as three separate fluids, loosely coupled. For 
the electrons and the ions, we assumed a perfect gas equation of state. 

Radiation Transport and Pressure 

Since the radiation pressure can range from being about equal to the matter 
pressure in front of the shock to more than 100 times the matter pressure behind 
the shock, it is important that this be treated properly. SIMPLE normally uses 
a full flux limited diffusion treatment of the radiation transport in the gas. For 
full diffusive radiation transport treatments to run in a reasonable amount of 
computer time, the zone size should be somewhere close to the actual optical 
mean free path. In a stellar interior, this mean free path is meters at most. With 
the large size of our area of interest our single zones were at least 1,000 km 
wide, so it was impossible to run simulations with full radiation transport in any 
available amount of computer time. Fortunately, this was not necessary. As we 
mentioned above, Weaver (1976) has shown that under the conditions present 
in supernova envelope shocks, the fluids should remain essentially in radiative 
equilibrium. Thus, we could mimic full radiation transport by assuming that 
the three fluids always maintained the same temperature and we modified the 
code to enforce this condition. To insure the optically thick approximation, we 
also set the opacity to a high value (lO^m" 1 ); a realistic, appropriate, calculated 
opacity would have varied over a large range, but the corresponding optical mean 
free path would have remained much smaller than our zone size, so we chose 
the fixed value for efficiency. This high value for the opacity helped to insure 
that the radiation remained in equilibrium with the other fluids. Our radiation 
treatment then consisted of an electron-radiation coupling routine, a flux-limited 
radiation pressure term, and a routine which summed the total energy present and 
apportioned it between the three fluids such that they all maintained the same 
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temperature. In order to check further both the equal temperature assumption 
and the general appropriateness of our overall radiation handling scheme, we ran 
the complete code on a finely zoned test case with full multi-group radiation 
transport. (Multi-group refers to the fact that we could break the radiation into 
frequency ranges which are treated separately by the radiation diffusion routine, 
to allow for departures from a purely Planckian local thermodynamic equilibrium 
distribution). For our test runs, the radiation was treated in 50 frequency groups. 
The zones were 100cm thick, approximately one optical Rosseland mean free 
path. The test was run in slab geometry for convenience. The entire test model 
was only 2,500cm thick. We input a shock with a velocity of 2.2 x 10 9 cm/s. 
Comparisons of the results for full transport and our approximation, at 6 x 10 7 5 
and 10 -6 s, are shown in figures 4 and 5 and the temperatures can be seen to 
agree to a very high degree (£ 3%). A detail of the three fluid temperatures near 
the shock front in the full radiation transport case is shown in figure 6. It can 
be seen that shock structure in the form of a small ion temperature spike only 
appears on a scale of ~ 10 - 20 cm (~ 1 - 2 zones immediately behind the shock 
front), very much smaller (by a factor of 10' 7 ) than our usual zone size. Further, 
the difference between the radiation and ion temperatures is only ~ 0.6%. Thus, 
we see that our equal temperature approximation is validated as predicted by 

Weaver (1976). 

Assumed Mechanisms of Importance 

The shock was mediated by an artificial quadratic viscosity term in the ion 
hvdrodynamics in the form of an added pressure Pq. 


Pq = f2.0p(Ax) 2 (gf) 2 , 

* Q l 0 


if zone is in compression, dV/dt < 0 
otherwise 


( 3 ) 


where p is the matter density, Ax is the zone size, V is the zone volume, and Au 
is the difference in velocity accross the zone. The multiplier of 2 is chosen by tra- 
dition. Test runs with other quadratic viscosity multipliers showed little variation 
in results apart from a slight rounding and smoothing effect at the leading shock 
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edge (see table 2). Full discussions of the Von Neumann-Richtmeyer artificial 
viscosity prescription in numerical hydrodynamics can be found in Larsen (19S9) 
or Mihalas and Mihalas (1984) (among others). The real physical mechanism of 
shock mediation in the regime of interest is radiation pressure, but ion viscosity 
was easier to implement in our code. The forced* equilibrium described above 
makes our method equivalent to radiative mediation (see also the discussion of 
analytic Hugoniot relations below for confirmation). Since our shock velocities 
are ^ 4 x 10 9 cm/s » 20 MeV /nucleon we are justified in using a non-relativistic 
^ ~ ' “ c ) assumption with respect to the overall velocities involved. More 

importantly, our post shock temperatures remain below 100 keV, so that our 
non-relativistic ideal gas equation of state remains valid. Another important ap- 
proximation is that we have assumed a lack of significant spontaneous electron 
positron pair production. Weaver (1976) has derived a criterion for whether pair 
creation processes are important ( i.e . when the number density of pairs equals 
the number density of ions): 


E 0 £ ,95$- 5 / 2 e* 
m e c 2 _ 5.94 x 10 9 

= ItT ” T\ 


( 4 ) 

( 5 ) 


where m e is the electron mass, is the immediate post shock temperature in 
Kelvin, and E 0 is the shock strength in MeV /nucleon. From this we can calculate 
that, since our shocks remain < 20 MeV j nucleon, our post shock temperatures 

must remain below 62 keV (which they do) if we are to be correct in neglecting 
pairs. 


We also included the effects of gravity in our simulations, but these turned 

out to have a negligible effect on the hydrodynamics and the nucleosynthetic 
yields. 
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Modeling the “Piston” 


To avoid the still unresolved question of the precise mechanism of shock 
production in core collapse supernovae, we placed a “piston” at the base of our 
envelope model and set it in motion with a test velocity. Its subsequent evolution 
was determined by the analytic models of Sedov (1959) for the propagation of 
a strong shock in a 7 = 4/3 gas (he. pure radiation) with a power law density 
profile. We changed power law exponents whenever the shock front entered a 
new region (figure 3) since the overall shape of the shock velocity distribution 
as it propagates will approach the shape characteristic of whatever the current 
power law is (Chevalier 1976). Thus, the shock front first- accelerates as it comes 
down the steep (w > 3) gradient of the edge of the helium zone, but then quickly 
(within ~ lOsec) becomes a decelerating shock as it settles into the u « 2 - 2.5 
region, which makes up the main body of the hydrogen envelope (and covers all 
of the narrow region of interest) (figures 2, 3, 7). We could also tell by inspection 
if there was an error in the velocity source since, after an initial energy transfer 
phase (~ 10-30 sec), the innermost zone (“piston”) should follow smoothly 
along behind the zones in front of it neither crowding them nor lagging them 
significantly. We could also observe any energy input or removed by an improper 
source function, after the initial energy transfer, by following the overall energy 
budget edits produced by the code (see below). 

In order to create a piston which will properly “follow along behind the 
post shock material, we used the Sedov similarity solutions for a propagating 
shock in a material with power law density falloff. As described above, we can 
approximate the density gradient in various regions of the star with a power law 
(see figure 3). The Sedov solutions then provide us with analytic formulas for 
the position of the shock front and a parametric set of equations from which we 
can derive the position of any zone behind the front at a particular time after 
the shock passes it. In practice, we characterize the strength of the shock by a 
parameter £ cx E 0 f ( 5 - u>) where E 0 is the total energy of a spherical shock blast 
in an infinite, perfect powerlaw density distribution of an ideal gas (initially 
unknown since the starting velocity is our independent variable), and w is the 
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density exponent as given in equation ( 2 ): 


£ — r 0 t Q . (6) 

Here to is the time it would have taken in a power law medium for the shock 
wave to travel from the blast center to the piston zone radius, r 0 . 

. _ 4 r 0 

° wo(w-5)(7+1)’ ^ 

where 7 is the perfect gas specific heat ratio. Since we start our simulation at 
tsim = 0 , the time, t in the subsequent Sedov equations is: 


t — tstm "b to (8) 

If, in addition, we let r 2 be the radial position of the shock front, v 2 be the (lab 
frame) velocity of the gas immediately behind the shock, and v be the velocity 

of the piston, then the Sedov equations for t and v in terms of the parameter V 
are: 


r 2 = r 0 g{ 7 ,w, V) 


(9) 


5— w 



i 


( 10 ) 


2 r 2 

(5 — w)f 


( 11 ) 


u = jW2(5-w)(7 + l)V/( 7 ,a;,V) 


( 12 ) 


where 
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(14) 
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(5 - A 7 

37 - 1 


2(6 - 37 - u) 
7(5 -w) 2 



_ 7 ~ 1 

012 ~ (u - 1)7 - 1 


(15) 


„ 7 = Ql l?7 A (16) 

6 (27 + 1 - 7 w) ' 6 - 37 - w 

Thus, we can now substitute V’s into these equations and get out pairs of 
(t si m,v) which we feed in as a rigid velocity source for the innermost zone of our 
mesh. We also keep track of the value of r 2 and, when r 2 (i.e. the position of the 
shock front) enters a new density powerlaw region, we take the t,v, and r at that 
point and use them as t 0 , v 0 , and r 0 for a new calculation with the u appropriate 
for that region. In these calculations, we have used 7 = 4/3 since the “radiation 
gas” is the actual primary mediator of the shock and the Sedov solutions are a 
“one-fluid” model and should thus be solved for the dominant “fluid . In oui case 
this is the photon gas (see Chevalier 1976). We should here note that model o58 
has an odd quirk in its structure with two power law regions with u significantly 
greater than 5 . In the other cases, such regions were small in terms of the time 
the shock spent in them and came at the start of the simulation only; thus, 
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approximating w with 4.9 posed no major difficulty. In this case, however, the 
presence of such large regions of high u> caused the code to be very unstable to 
“overdriving” errors, so we attempted to get a reasonable piston velocity profile 
by “turning off” the piston (he. setting the velocity to a value of about 1/10 of 
the starting velocity) after 50 seconds. The shock continued to propagate nicely, 
but with some pinning down’ of the shock tails as in the linear piston case 
(Figure 12). Since the post-shock temperature was still well below the critical 
temperature of ~ 10 keV this should not change the final yields by a great deal. 

The Reaction Network 

The network of reactions we used is shown in table 3. The thermonuclear 
reaction rate values were taken from Coughlan and Fowler (1988). We chose 
this particular set of reactions by examining their tables of rate versus temper- 
ature. We began with the fundamental reaction of interest, 3 jfe(a, j) 7 Be, and 
then followed each possible reaction with the ambient nuclides (a, 3 He , H) until 
each possible branch hit a “wall”, i.e. a reaction which would not occur below 
~ 100 keV . We then cast this network into an implicitly differenced routine which 
would evolve the network with time. We found that at our temperatures of inter- 
est, the network was stable up to timesteps of 0.5s. To justify the completeness 
of picking just these reactions, we used our implicitly differenced mesh to recom- 
pute the generic explosion models of Arnould and Norgaard (1975). They used a 
complete light element network containing 42 nuclei linked by 178 reactions and 
P decays. For a model temperature vs time profile they assumed Arnett’s (1969) 
generic explosive nucleosynthesis formulation, in which one assumes that an ideal 
gas is instantaneously compressed to an initial density, p l and temperature T 1 ,. 
The gas then exands adiabatically so that: 

P = Pit ~ t,T ” (17) 

T = T t e~ </3rei (18) 

where r ex is a characteristic gravitational free expansion time scale: 
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(19) 


Tex = X 446 Po 1/2 

and x is considered a free parameter. Figure 8 shows the. results of running our 
nucleosynthesis network using the same conditions as Arnould and Norgaard. A 
comparison with their figure 4 reveals only small differences, with the exception 
of the final 3 He abundance which is a factor of ten lower. This is due to the 
higher 3 He (0,7 ) 7 Be rate used in our network, which leads to a more rapid de- 
pletion of 3 He in our case and consequently a lower final abundance. This rate 
was recently revised as the result of new experiments (see Coughlan and Fowler 
(1988) and Kawano, Schramm, and Steigman (1988) and references therein). We 
also compared our numerical mesh (without the nonlinear 3 He+ 3 He term) to an 
analytical solution of the differential equations (solved for a fixed temperature) 
obtained by the Laplace transform method. We used the same starting compo- 
sition of nuclides as in the Arnould and Norgaard simulations. Figure 9 shows 
the two models and, except for the slight irregularity in the destruction curve 
of 11 B at times of ~ 10 -5 to 10 -2 s, they are identical. As a point of interest, 
inclusion of the 3 He + 3 He reaction results in a slightly faster destruction of the 
3 He starting material and corresponding lowering of the yields of all of the othei 
products (see figure 10). 

A close examination of the results of these simulation will disclose the overall 
path of C ameron- Fowler lithium synthesis. Figure 10 illustrates the reactions 
involved. The basic paths are: 


He+ 4 He S Be + 7 

(20) 

Be +* He — > n C + 7 

(21) 


The 7 Be and n C will become 1 Li and n B by electron capture and positron 
emission, respectively, as the ejected material cools. Any initial 11 B is quickly 
destroyed by protons. Interestingly, this whole scheme would be phenomenally 
unproductive if it weren’t for the photodisintegrations. One would expect that 
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proton destruction of both 7 Be and 11 C into *B and 12 N should be disastrous for 
this mechanism, since both S B and 12 N very quickly /3 decay into 2 a’s and 12 C 
respectively. However, at temperatures above 1.2 x 10 8 for 8 B and 4.8 x 10 8 for 
12 N the photodisintegration of these elements back into 7 Be and 11 C takes place 
at a faster rate than beta decay. This behavior is evident in figure 8, where first 
the 12 N and then the S B show fairly sharp “bumps” as the temperature falls off. 
This is due to a rapid rise in their abundances as the photodisintegration becomes 
negligible followed by a decrease in abundance as the production reactions for 8 B 
and 12 N also cut out with falling temperature, thus allowing them to beta decay 
away. This same basic pattern of events will characterize any Cameron-Fowler 
type light element synthesis environment: that is, any one characterized by a 
sharp, large temperature increase followed by rapid cooling. 


Quality Control 

After glancing over the Sedov equations for a propagating shock, the first 
quality control question the reader may be tempted to ask is just how sensi- 
tive the light element production is to the piston velocity curve. To test this 
qualitatively, we ran three simulations of approximately the same final energy 
(as determined by the energy audit of the simulation code) but with radically 
different piston velocity curves (shown in figure 11). The “correct” curve was 
formed by the method described above. The “bad” curve was formed by the 
above method as well, but a much higher starting velocity was combined with 
the wrong density power law exponents, a/s (we leave it to the reader’s imagina- 
tion to discern the circumstances of this simulation’s generation). The “linear” 
piston function was produced merely by drawing a line between two of the points 
of the correct curve. Velocity and temperature profiles for the three pistons are 
shown in figures 12 and 13, and except for the “pinning down” of the piston 
velocity in the linear case resulting from “overslowing” the piston (an opposite 
effect occurs if the gas is “overdriven” by a piston moving too fast: velocity and 
temperature “tails” develop) , there is not a great deal of difference between 
them. As can be seen from table 4, the 7 Li yields are similar (within 34%) and 
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even the amount of kinetic energy error (see below) is similar. The difference 
between the 11 B yield of the “bad” piston and the other two is a result of the 
very high initial velocity, which causes a relatively large build up of 11 B in the 
zones immediately adjacent to the piston zone because of the lack of hydrogen, 
i.e. destructive protons, in that region. This production site is not very signifi- 
cant if there is even reasonably efficient production of 11 B in the same zones as 
the 7 Li production, since 3 ife is already being turned into 7 Li and 11 B at 100% 
efficiency in this region. Thus, more favorable shocks can increase 11 B production 
farther out, but the production at the edge of the 4 tfe zone remains fixed. For 
example, in the case where we ran a trial version of the code with zero radiation 
pressure resulting in very high post-shock temperatures, the contribution of the 
inner zones to the total 33 B production was small (less than 5%). 

In terms of overall quality control, we used two methods to determine whether 
the shock structure was accurately modeled. First, we compared our shock struc- 
tures at given instants to the the general effective photon shock model Hugoniot 
relations of Weaver (1976), which describe the “jump” in the values of of the 
hydrodynamic variables across a shock front of arbitrary strength. In Weaver s 
non-dimensional form these are: 


X<i r lf T f + 2 Tf = + e 0 + 1) 


( 22 ) 


3 X?7/ T / + r /( a e + 3 / 2 ) = 1 / 2 ( 1 _ w)( 1 - Vf + 2e °) + a ° e ° ( 23 ) 

where: 
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In terms of our parameters these can be rewritten in the form of the continuity 
equation 


u gas 

K = 1 - — 

U o 

the momentum conservation equation 


Po 

P l 


(27) 


2kR P0 T 1 + — T 4 = p 0 u 2 0 (l - 1/ k) + T { 4 + 2Rp 0 T 0 , (28) 

and the energy conservation equation 


(2a e + 3 )k.RpqTi + 6— T 4 = po«o«(l ~ l/«) 2 + 2 |^( 4/c - l)?o 4 + 2 P0 RT 0 (5k - 2). 

(29) 

oc e = 3/2( ^ 373xl0l 8 + 1) (30) 

T f + 5 

In these expressions, To and po and T\ and p\ are the temperature and density 
immediately in front of and behind the shock, respectively, no is the speed of the 
shock front relative to the unshocked gas, and iig as is the (lab frame) velocity 
of the gas behind the shock. All quantities are in cgs units and temperature is 
in K . Readers familiar with the totally opaque gas, strongly radiative shock 
approximation (see t.Q. hlihalas and hlihalas (1984)) will note that this energy 
equation is different. This is due to the inclusion of terms accounting for the 
energy flux across the shock front. Thus, while a completely ion mediated shock 
(i.e. no radiation pressure) will produce spuriously high temperatures behind the 
shock, what one might naively consider a fully radiating shock creates post-shock 
temperatures which are a factor of ten or so too small. Table 5 shows the compar- 
ison of our simulation and the Hugoniot relations for a wide variety of conditions 
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which span those encountered in our models. The error is generally between 5 
and 10 % with a few larger errors appearing non-systematically. The greatest 
source of error here (easily up to 10% in any case) is due to the impossibility 
of actually locating the shock front in an artificial viscosity code. To measure 
the speed of the shock front (our shock strength defining parameter), we noted 
every few cycles of the code which zone was controlling the time step. This zone 
corresponded to the one with the greatest compression and so “contained the 
shock front at that time. By dividing the original width of the zone by the time 
the shock spent in it, we obtained an estimate of the speed of the shock. Because 
of the finite size of these zones (each was 100 cm initially), the finite size of the 
timestep, and the spreading of the front over more than one zone in the weak 
shock cases, this could lead to significant “error” which does not reflect an inac- 
curacy of the actual hydrodynamics of the code. The other source of error was 
closely related. Since the zones behind the shock rang (see above) and because 
the shock is spread over a small number of zones it was often difficult to deter- 
mine which zone to use for the actual peak post-shock velocity and temperature. 
In some cases where we checked, there was as much as 10 % difference between 
adjacent zones behind the shock. Considering these uncertainties, the results in 
table 5 indicate good agreement between our code and the analytic predictions 
of the Hugoniot relations. 

The second check was to track the overall energy flow between the thermal 
energies of the three fluids (electrons, ions, and radiation) and the kinetic energy 
of the shocked material. By comparing the total kinetic energy created by the 
source plus the initial thermal energy of the mesh with the total energy present 
in the mesh at any time, we could tell if there was any significant unaccountable 
loss or gain in energy. The generally accepted rule for artificial viscosity codes 
is that overall energy error should not exceed ~ 10 - 20%. By far the largest 
errors occur in kinetic energy, while the internal energies of the electrons and 
ions always remain much smaller than the kinetic energy error. A large part of 
this “error” is created by the way the hydrodynamic equations are transformed 
into difference equations. For stability reasons, the algorithm is structured so 
that the position of a zone is correct at the end of the time step, but the velocity 
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is off by approximately one-half of a time step. Thus, when the code calculates 
the current kinetic energy, it uses the velocity at the previous half time-step. If 
the velocities are changing rapidly (which is true at the shock front) this can 
create an error which looks like missing kinetic energy. Our errors ranged from 
~ 20% to ~ 4% with few exceptions. These were all in the very low energy 
shocks in the heavy blue stars. This error comes from at least two sources. First, 
the post shock ringing discussed above is much closer in magnitude to the actual 
shock size in these very weak shock situations. Further, gravitation (which was 
added on to the finished code, and thus not included in the energy conservation 
routines) while having an almost imperceptible effect (~ 1%) on light element 
yields, did cause a comparatively large “error” in the kinetic energy of these 
weak shocks. It should also be noted that these errors build up with time since 
the error per step remains fairly constant (with some decrease as the shock gets 
weaker), and are much smaller at the end of nuclear burning around 100 seconds 
than they are at the very late time of 1,000 seconds reported for most of the 
simulations. For instance, model m251 with shock energy 3.05 x 10 50 er^s at 150 
seconds has deposited all of its piston energy in the outer envelope (i. e.the piston 
is coasting’) and has finished both 11 B and 7 Li production, but it has only a 
14% kinetic energy error as opposed to 38% at 1,000 seconds. 

As a final qualitative check we “kept an eye on” the overall shape of the shock 
as it evolved to make sure no “tails” or “pinning down” occurred, as in the linear 
piston case shown in figures 12 and 13 and in the o58 model as discussed above. 


Results and Conclusio ns 

Our overall Li and yields for a range of energies and precursor masses 
are shown in figures 14— 19and in table 6. Most of the simulations were run for 
1, 000 simulation seconds at which point the maximum shock temperature was 
well below the minimum temperature for nuclear processing of 1 x 10 7 0 K. The 
exceptions were the very high energy shocks run on the m251 model and the 
the very massive models (o58 and olOO), which were run only until the shock 
left the region of interest (z. e.the 7 Be and U B abundances stopped evolving) 
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because of cpu time constraints. The time evolution for a typical shock passing 
through the region of interest is shown in figure 20. The 1 Li and 11 B were 
produced exclusively in a region centered about the base of the surface convection 
zone (within the first solar mass worth of material outside the He coie shell) 
where the hydrogen and helium have approximately equal mass fractions (see 
figures 21 and 22). For purposes of comparison, in other detailed explosion 
simulations and in observations the total energy deposited by the shock is found 
to be 1 - 3 x 10 51 ergs. It should be born in mind that the energy we report is 
deposited in the outer envelope and does not include energy residing in the ejected 
core. In table 6, we also give the total kinetic energy present at the stopping 
points for our simulations. We note that at these points in time, the shock is 
still deep within the star and therefore most of the shock energy still resides in 
the radiation immediately behind the shock; this implies that the energy of the 
core ejecta at the end point of the simulation is comparatively small and thus our 
reported “total energy” is a good number to compare to the measured kinetic 
energy at light breakthrough in observations. The shock weakens as it progresses 
outward to times past the end of our simulations, and thus a progressively smaller 
percentage of its energy resides in radiation; by the time the shock reaches the 
surface of the star, most of its energy is thus kinetic (Chevalier, 1976). 

As expected, the red giant stars (models ml51,ml52) are much poorer pro- _ 
ducers of light elements than are the blue giant stars (all of the other models). 
The exception to this is the massive star model with mass loss which is a very 
poor 7 Li producer because it lost most of its 3 He as a consequence of mass loss, 
and thus had a much lower preexplosion 'Hit abundance than the other models 
(see table 1). It is also clear that only at very high shock energies, ^ 4x 10 51 ergs , 
in our most favorable models involving large blue giant pre-supernova stars, does 
this process even reach a Pop I 1 Li abundance overall in its ejecta. For energies 
(~ 10 times) higher than observed, the 1 Li mass fraction averaged over the total 
ejected mass of the star is still only 5 x 10 9 . 

Figures IS and 19 show the total mass of 1 Li and 11 B produced as a func- 
tion of stellar mass for the most productive pre-supernova models (including the 
unrealistic olOO). There is not much variation in production over this wide mass 
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ran § e - Thus, it seems unlikely that there is a hidden w super- producer” class of 
stars. We shall discuss below how to recognize such a class of models if they do 
appear in future research, however. 

As derived in DSST, we may compare our yields to the expected iron pro- 
duction of core collapse supernovae. As argued by Arnett, Schramm, and Truran 
(1989), because the central cores of all Type II supernovae are so similar, we 
might expect them all to yield comparable quantities of iron to that determined 
for supernova 1987A, or ~ 0.07 Mq of iron. Type I supemovae can also account 
for a significant fraction of the galactic iron. Thus, to be a useful source of 
Li and 11 B, the ratios Li/F e and BjFe in the ejecta must be greater than or 
comparable to present values in the Galaxy. DSST then use these numbers to 
obtain the result that, to within factors of a few, a w critical” yield of lithium in 
a supernova is 


log M(Li ) = -6.3 ± 0.1 M(Li) <7 4 x 10 -7 Mq. (31) 

Similarly, for boron they find, 

log M(B ) = -6.5 ± 0.1, M(B ) £ 3 x 10 -7 Mq. (32) 

However, the lithium vs. iron measurements also indicate that over a range 
of iron abundance (mass fraction) (Fe/H) from 10~ 4 to 10" 1 - 5 times solar, the 
lithium abundance does not increase by more than a factor of 2 (i.e.an increase 
of ~ 10 in (Li/H)). This implies that the average lithium output ratio from 
any Type II supernova process must obey 

Li IQ — 10 r * 

(pT at ~ 1q9 ^) ^ 2Q— 1.5 ( jjTj now ) ~ 10°' 5 10 -6 - 35 10 -515 . (32) 

If we further assume that these supernovae produce ^ . 1 Mq of iron, we obtain 

M(Li) £ 10 -5 ' 15 10 -1 - 7 x 1O _7 M 0 . (32) 

Even if we allow for the very approximate nature of the assumptions in these 
calculations, they indicate that effective supernova production is rather tightly 
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constrained to a range of about 10 7 to 10 ® Mq of 7 Li per event averaged over 
all masses. This constrained range argues for a slightly “slower” production site 
such as high mass AGB red giants (see below). 

As can be seen in our models, pushed as far as we think reasonable (and 
about as far as we can trust our code), 7 Li yields are still a factor of 80 below 
this “break even” point and the 11 B yield is orders of magnitude too low. 

It, thus, seems unlikely that the envelope shock mechanism can independently 
account for the evolution of 7 Li abundances from Pop II to Pop I values. It is 
important, however, to try to ascertain how firm this conclusion is. In terms 
of model sensitivity, we have two questions to answer. First, given a potential 
production site (t.e.gas containing 3 He ), what temperatures and densities are 
most likely to produce significant amounts of 'Li and 11 B ? Arnould and Nor- 
gaard (1975) have investigated this in detail, so we can use their results as a 
guide to determine when more detailed modeling might discover a meaningful 
production site. Second, since their results require the input parameters pi and 
Ti (equations (17)- (19)) we must ask what explosive situation will produce this 
necessary density and temperature. In terms of our envelope shock model, we 
can use Weaver’s (1976) Hugoniot relations for hydrogen envelope conditions to 
tell us what input shock conditions will produce significant 'Li . Figure 23 de- 
tails the range of interest for our type of scenario. Post-shock conditions do not 
depend strongly on pre-shock temperature. We chose 6 keV as a reasonable value 
for all of the possibilities, since at temperatures much above this no 3 He would 
survive anyway on stellar timescales. Then we fed in a range of shock velocities 
and starting densities. The different symbol types in figure 23 represent different 
shock speeds. Each point in figure 23 represents an increment of 1 gm/cm 3 over 
its neighbor of the same symbol type. The area to the upper right of the line, 
’limit’, represents Arnould and Norgaard’s region of efficient 7 Li production for 
their most productive model. The greatest sensitivity (in terms of conditions 
that might really be met in stars) is to density variations. At 10,000 km/s, 
a maximum reasonable shock speed value in current models, the threshold for 
significant ‘Li production occurs around p = 7 gm/cm 3 . This condition is def- 
initely not met in any current envelope models that we know of. We conclude, 
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then, that any likely site for shock-induced production of 7 Li would require the 
presence of 3 if e with a mass fraction of ^ 10 -7 and either a preshock density of 
7 g/cm i 3 or a shock strength 20, 000 km/s ~ 1.8 MeV /nucleon. 

In a subsequent paper (Brown 1990), other non-cosmological sources for 7 Li 
will be investigated in detail. It may nevertheless be useful to mention briefly 
the two alternative non- super novae mechanisms for galactic 1 Li nucleosynthesis 
which are considered to be promising candidates: those involving classical novae 
and red giant stars. 

Classical novae can produce 1 Li . The critical question is whether they 
can form sufficient quantities to perhaps make an important contribution to the 
abundance of 7 Li in the Galaxy. The mode of 7 Li production here is similar 
that in red giants: 7 Be is carried outward by convection to cooler regions of the 
envelope on a sufficiently rapid timescale to ensure that it will not be destroyed 
via 7 Be (p, j) 8 B . The total mass of 7 Li ejected is a sensitive function of the 
conditions achieved in the outburst, particularly the convective history of the 
ejected matter, and may thus be expected to vary somewhat from event to event. 
The main problem is the fact that the 7 Li overabundance factors are less than 
those of other nucleosynthesis products. The hydrodynamic models of Starrfield, 
et. al. (1988) predict average enrichment factors for the rare isotopes 13 C and 
15 N which are approximately 10 times that of 7 Li . If these relative production 
ratios are correct, and are representative of classical nova systems, it must follow 
that novae can account at best for only 10 % of the 7 Li in the Galaxy - when 
one assumes that 13 C and 15 IV have their origin entirely in nova explosions. 

The observation of high ' Li concentrations in luminous red giant stars (Smith 
and Lambert, 1989) suggests that these may also play an important role in the 
production of 7 Li in the Galaxy. The difficulty here is to obtain a realistic theo- 
retical estimate on the basis of our current imperfect knowledge of the evolution 
of the surface compositions of asymptotic giant branch stars. The possibility 
that red giants may be significant contributors can, however, be demonstrated 
in a model independent manner. Following the discussion of Scalo (1975), we 
assume that lithium-rich red giants constitute a fraction 1/100 of the “peculiar” 
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red giants and that these peculiar giants, in turn, represent a fraction 1/100 of all 
red giants. Since red giants represent approximately l/10th of all stars, we then 
estimate that there are of order 10 10 red giants. The bright, luminous red giants 
that show enhanced lithium concentrations are estimated to be losing mass at 
rates approaching 10 — 7 solar masses per year; such stars also show lithium en- 
richments of a factor of order 10 to 30 relative to solar system matter. If we take 
these factors together, we arrive at the following rate of production of lithium in 
the Galaxy: 

— — 10 10 10 -7 1CT 8 = 10 ~ 9 M Q /yr (33) 

100 100 

If we now assume that this continues over 10 10 years of galactic evolution and 
contaminates approximately 10 10 solar masses of gas in the ISM, we see that red 
giants can roughly account for the observed 7 Li abundance of 10 9 by mass. It is 
clear that a detailed calculation is required that takes properly into account the 
nucleosynthesis conditions that are believed to obtain in red giants. However, 
the above argument indicates that we must at least be aware of the possibility 
that red giants can synthesize lithium. An additional factor to be considered is 
the fact that they must be able to do this on a relatively rapid timescale, early 
in galactic history. 

In conclusion, Steigman (1990) has done some analytic calculations of lithium 
evolution in the instantaneous recycling approximation and the preliminary con- 
clusions seem to be that if the lithium “tracks” the iron abundance, it fits the 
observed data better than if it tracks the oxygen abundance. This would seem 
to imply that high mass rather than intermediate mass stars are the most likely 
sites for lithium production. If this is born out by more detailed models, it could 
point to the need for an effective mechanism of 7 Li production in supernovae if 
the Population II lithium abundance is primordial. 
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TABLE CAPTIONS 


1. Model parameters include total mass, envelope mass, hydrogen in envelope, 
final convection zone 3 He abundance, and the mixing length used in stellar 
evolution. 

2. Results and yields for different quadratic artificial viscosities. These are 
all model m251 with the same initial input velocity and audited at 200 
simulation seconds. There is very little difference between the simulations. 
Graphs of the hydrodynamic variables show a slight softening (rounding 
over) of the leading edge of the shock front in the q — 5 case. 

3. Table of the reactions included in our thermonuclear reaction network. 

4. Energetics and final abundances of 11 B and 7 Li for the three pistons. 
Notice that there is scant difference in results except for the higher 11 B 
abundance in the flawed piston model which is discussed in the text. 

5. This table shows our systematic check of the numerical code against the 
supernova envelope shock Hugoniot relations of Weaver (1976). The percent 
error is rarely more than 10% and only once over 20%. These errors are 
all attributable to the uncertainty involved in locating the actual shock 
front and deriving its speed as discussed in the text. The columns give 
the pre-shock density and temperature (To,/>o), the velocity of the shock 
front itself {v s h oc k), and the post shock velocity and temperature (vi,Ti) 
calculated analytically with the Hugoniot relations and evolved with the 
numerical code. 

6. This table presents our overall results. The columns give: mass, source 
energy, kinetic energy, kinetic energy error, simulation seconds run (since 
errors grow with time), ^ Be yield and 33 B yield. 
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FIGURE CAPTIONS 


1 . This well-known big bang nucleosynthesis graph shows the primordial abun- 
dances of the big bang produced elements vs. the average density (today) of 
the Universe. The “allowed window” is derived from present limits on the 
abundances of these elements. Note how the observed Pop II abundance 
sits right in the dip between the direct 7 Li regime and the 1 Be regime. 

2. A schematic cross section of the pre-supernova star. The dotted line, marks 
the lower edge of the surface convection zone which is the main area of 
lithium production. The drawing is to scale for our model m251 for all but 
the surface of the star which should be roughly 3 times farther from the 
center. Radius represents mass enclosed. 

3. a) Plot of log(p) vs log (radius) for m251 with the various power law regions 
marked. We did not compute u/s for more than three regions since our 
simulations had run their course (i.e.the post shock temperature was much 
too low for nuclear processing) well before the time the shock reached the 
outer edge of the third region. The slopes and their statistical accuracies 
are: region I- u> = 5.5 ± .4, p — .97; II- u = 2.3 ± .06, p = .99; III- 
<*; = 1.713 ± .006, p = .9997. Since the Sedov (1959) solution is only valid 
for u < 5 we approximate u for the first, relatively small, region by u = 4.9. 
b) A plot like “a” shown for model o58, note the two areas of high (> 5) 

u. 

4. Comparison of ion temperature profile with full radiation transport and 
our approximation using enforced equilibrium of the three fluids at 6 x 
10 -7 seconds and 1 x 10~ 6 seconds. 

5. Comparison of velocity profile with full radiation transport and our ap- 
proximation using enforced equilibrium of the three fluids at 6 x 10 _7 s and 
1 x 10 -6 s. 

6. Closeup of the three fluid temperatures vs. radius in the radiation transport 
test case at 6 x 10~ 7 s. Note the essential lack of deviation from complete 
equilibrium. 
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7. Velocity vs. time for the piston for model m251, starting velocity 12.6 x 
10 8 cm/s (final energy = 1.24 x 10 51 ergs). Note the rapid initial increase 
in piston velocity as the shock “rolls down” the gradient at the edge of the 
4 i7e zone, followed by a sharp decline in speed as it begins to sweep up 
the material in the “flatter” (i.e.lower density gradient) region of the star. 

8. This log-log plot of abundance vs. time for all of the active nuclides in our 
network using a x = 1, T 0 = 4.5 x 10 8 °K, po = .911 generic explosive nu- 
cleosynthesis model can be compared with Amould and Norgaard’s (1975) 
figure 4. 

9. Plot of mass fraction vs. time for all of the active (i.e.evolving) nuclides 
in our mesh for the single temperature analytic solutions and the numeri- 
cal mesh r both computed without the '’‘He A 3 He reaction. Temperature 
equals 3 x 10 8 °K. The starting abundances are those chosen by Amould 
and Norgaard (1975). 

10. Plot of the full numeric solution with the 3 He A 3 He reaction for the con- 
ditions of figure 9. 

11. Velocity vs time for correct tu, flawed tu, and linear pistons. The flawed lo 
piston starts at velocity 3.78 x 10 9 cm/s. 

12. Velocity vs radius at 50 second intervals for the three pistons shown in figure 
11. Note how the slow speed of the linear piston in the later snapshots “pins 
down” the tail of the velocity making the linear piston velocity profiles look 
obviously different from those of the other two pistons. 

13. Temperature vs radius at 50 second intervals for the three pistons. 

14. The overall output mass fraction (roughly) of 7 Li is given as a function 
of total input shock energy. The mass fractions are found by dividing the 
total mass of 7 Li produced in a model (see table 6) by the total starting 
mass of the star minus 1.4Mq. 

15. This is the same as figure 14 for 11 B production. 

16. This shows the same as figure 14 but only for model m251 and includes the 
very high energy shock data. 
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17. This is the same as figure 16 for 11 B . 

18. This shows the total mass of 7 Li produced as a function of initial star mass 
for various total shock energies. The models chosen for each mass were those 
which produced the most 7 Li (this includes the naive model olOO with no 
mass loss rather than the probably more realistic o58). 

19. This is the same as figure 18 for n B . 

20. Temperature vs. radius and velocity vs. radius for the main nucleosynthesis 
part of the simulation for model m251, energy 1.24 x 10 51 ergs. The curves 
represent successive “snapshots” at 10, 20, 30, 40 and 50 seconds from start 
of simulation. At 50 seconds, 99.3% of the final energy of the shock has 
been transferred from the piston to the mesh, 80% of the final 7 Be has been 
produced, and all of the final 11 B has been created. 

21. Graph of 7 Li mass fraction after shock processing vs mass enclosed (mass 
shell) for model m251, shock energy = 1.24 x 10 51 ergs. Note that all of 
the significant production occurs in a region in mass shell coordinates of 

1 M© in width. 

22. This is the same as figure 21for 11 B . As described in the text, the poor 
conditions for 11 B production in these shocks cause the 4 He shell edge pro- 
duction to greatly overshadow the “main production region” which would 
dominate under more favorable temperature-time profiles. The small fea- 
ture at the left edge of figure 21 (at ~ 7 Mq) corresponds to the extremely 
large feature at the left edge of this figure. 

23. This graph plots log(pi) vs T,/10 8 for different shock strengths. The line, 
’limit’ represents the edge of the region of significant lithium production 
in Amould and Norgaard’s (1975) most favorable scenario (i.e.x = 1-0 in 
our equation (19)). Each symbol type represents a different shock strength 
(the number after “v” refers to the speed of the shock front in thousands 
of km/s). The lower left point in each “v” curve represents a pre-shock 
density of 1 gm/cm 3 and each point to the right of that represents an 
increment of 1 gm/cm 3 . Thus, the curve vlO crosses the limit line on the 
seventh dot meaning that a shock with velocity 10, 000 km/s will efficiently 
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bum any *He present into 7 Be if the gas has a density 7 gmjcnr? before 
the shock hits. Since the post shock conditions are fairly insensitive to pre- 
shock temperature, all of the points are calculated for a single pre-shock 
temperature of 6 keV . 
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TABLE I 


Model Name 

Mass 

(Solar Masses) 

Envelope Mass 
(Solar Masses) 

Hydrogen Fraction 
of total envelope 

X(He3) 

Maximum" 


ml50 

15 

8.951 

.648 

1.032E-4 

1.4 

ml51 

15 

9.647 

.625 

9.450 E-5 

EQ 

ml52 

15 

9.624 

.622 

9.441E-5 

1.6 

m251 

25 

14.995 

.570 

8.893E-5 

1.4 

m252 

25 

15.009 

.555 

8.893E-5 

EQ 

o58 

58. 6 C 

11.3 

.170 

6.077E-6 

1.6 

olOO 

100 

54.7 

.435 

6.373E-5 

1.6 


a The initial input 3 He abundance, 2.1 x 10~ 4 , is the same for all models but convection 
has depleted this abundance. The value in this column is the amount left at the surface 
of the star (and is approximately constant throughout the convective region of the 
envelope) immediately before explosion. 

6 a is the ratio of the mixing length to pressure scale height. 

c This is the final preexplosion mass of the star. It began life as a 100 star but 
suffered mass loss.- 


TABLE 2 


Viscosity 

multiplier 

source energy | kinetic energy 
(10 51 erjs) 

kinetic energy 
(% error) 

7 Li | 11 B 

(solar masses) 

0.5 


1.78566 

3.9779 

1.58940e-09 

9.92745e-16 

1.0 


1.80532 

3.0472 

1.32150e-09 

7.90642e-16 

2.0 

4.08258 


2.4725 

1.11620e-09 

6.29119e-16 

3.0 

4.06808 

1.86078 

2.2276 

9.87200e-10 

5.20774e- 16 

5.0 

4.04182 

1.88856 

1.9804 

8.17190e-10 

3.81102e-16 


TABLE 3 


3 He + 3 He — > 

0“ 


3 He + a — * 

7 + 7 Be 

— (3 + 7 Li h 

7 Be + p < — ► 

7 + 8 B 

— ^ + 0 

7 Be + a — ► 

7 + 11 C 

~^f3+ n B 

n C + p «— » 

7 + 12 N 

— /? + 0 

n B + p — . 

0 


n C + a — +■ 

0 



a A 0 indicates that the reaction products are not accounted for in our reaction network. 
E. g. any production of 4 He would be insignificant given the enormous amounts of it already 
present, so any reaction which produces only 4 He enters only as a destruction reaction. 

b j3 refers to any nuclear process involving electrons. The electron capture rate for 7 Be + 
e~ — y 7 Li is so slow r x / 2 ~ 49 days that we set it to zero in the actual equations. 


















































TABLE 4 


Piston 

source energy 

kinetic energy 

kinetic energy 


rr B 

name 

( 10 51 er< 7 $) 

(% error) 

(solar masses) 

line 

1.50076 

0.45983 

8.5942 

3.8427e-10 

3.0647e-18 

correct 

1.44345 

0.41938 

9.4794 

3.6502e-10 

3.2578e-18 

flawed 

1.74107 

0.32704 

12.668 

2.5151e-10 

3.4536e-16 


TABLE 5 


To 

(keV) 

0.1 


P o 


(g/cm 3 ) 


0.01 


5.0 


^shock 

(10 8 cm/s ) 

L285 

6.185 

11.98 

23.02 

1.466 

6.341 

12.42 

23.68 

1.371 

6.330 

12.57 

26.19 

2.244 

6.640 

12.85 

28.54 

36.40 


analytical 

1.012 

5.255 

10.24 
19.70 

1.007 
5.333 
10.58 

20.25 
8.880 
5.275 
10.67 

22.38 
1.103 
5.330 
10.75 
24.32 
31.08 


(10 s cm/s ) 
numerical 
1.005 
4.969 
9.408 

18.39 
1.013 
5.074 
9.935 

19.22 

10.07 

5.089 
10.10 
19.67 

1.022 
5.105 
10.12 
20.02 
28.88 


% difference 
(L69 
5.59 
8.46 
6.87 
0.59 
4.97 
6.28 
5.21 
12.5 
3.58 
5.48 
12.8 
7.62 
4.31 

6.03 

19.3 
7.33 


analytical 

0.992 

2.845 

4.039 
5.640 
1.882 
5.021 
7.256 
10.14 
2.352 
8.569 
12.80 

18.89 
6.403 

12.89 

19.07 

29.39 
33.34 


Ti 

(, keV ) 
numerical 
0.968 
2.650 
3.110 
5.477 

I. 890 
4.611 
6.195 
9.792 
2.124 
7.725 

II. 27 
17.64 
6.515 
13.94 
18.01 
26.97 
32.47 


% difference 
2A4 
tT09 
25.9 
2.93 
0.42 


10.1 

10.3 

12.7 

6.84 


2.64 






































TABLE 6 


Model 

name 

source energy | kinetic energy 
(10 51 ergs) 

kinetic energy 
(% error) 

simulation 

seconds 0 



ml50 

0.29385 

0.15486 

15.855 

1000 

3.2158e-12 

1.2333e-22 

0.97105 

0.37325 

10.219 

1000 

1.9340e-ll 

7.3619e-19 

2.18326 

0.72148 

9.2790 

1000 

5.6805e-ll 

4.2014e-17 

3.87165 

1.12993 

9.2871 

1000 

1.1455e-10 

4.7348e-16 

ml51 

0.34535 

0.31627 

6.9998 

1000 

1.6752e-14 

9.2134e-29 

1.38906 

1.14966 

6.4274 

1000 

4.5096e-13 

1.6257e-24 

2.89712 

2.25534 

6.4601 

1000 

1.8693e-12 

1.2783e-22 

3.53510 

'2.65417 

4.6897 ' 

1000 

3.2910e-12 

1.3568e-21 

ml52 

0.33863 

0.30640 

7.5224 

1000 

8.5421e-14 

6.3490e-28 

1.32960 

1.08162 

6.9534 

1000 

2.1383e-12 

8.4426e-24 

3.02116 

2.30563 

6.4905 

1000 

8.6564e-12 

6.1905e-22 

3.75584 

2.69552 

4.9383 

1000 

1.5630e-ll 

6.6912e-21 

m251 

0.30519 

0.10868 

38.051 

1000 

2.0853e-ll 

1.0467e-22 

1.24211 

0.28112 

20.664 

1000 

2.5416e-10 

7.2309e-19 

2.76511 

0.50675 

' 15.343 

1000 

8.9687e-10 

6.5810e-17 

4.13690 

0.63388 

14.100 

1000 

1.7477e-09 

6.2910e-16 

6.64756 

1.50023 

6.051 

147 

3.3559e-09 

3.7072e-15 

9.12673 

1.95884 

5.715 

130 

4.9402e-09 

1.3272e-14 



0.10225 

38.803 

1000 

1.2461e-ll 

6.8521e-25 j 


0.25737 

20.814 

1000 

1.5904e-10 

5.8135e-22 

2.20954 

0.46378 

15.024 

1000 

5.8167e-10 

7.3852e-20 


0.60426 

13.313 

1000 

1.1651e-09 

2.2900e-18 

olOO 

0.61285 

0.50398 

19.1328 

200 

8.2099e-12 

6.4412e-27 

2.81542 

1.99089 

6.2674 

100 

1.0315e-10 

8.4875e-24 

5.57902 

3.44104 

5.6527 

150 

3.1250e-10 

2.8680e-22 

9.74379 

5.62623 

4.9346 

150 

6.6917e-10 

2.8020e-21 

o58 

0.93739 

0.18375 

179 

700 

1.4003e-13 

1.08214e-27 

3.10352 

1.87221 

28.2923 

650 

1.8622e-12 

9.9448e-26 

16.0503 

7.49521 

7.4246 

350 

3.9765e-ll 

4.6342e-23 


° Energetics errors tend to grow with time. 
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